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ABSTRACT 

Presents  a  procedure  for  computation  of  yield  tables  for  diseased, 
managed,  even-aged  stands  of  lodgepole  pine  in  Colorado  and  southern 
Wyoming.  Stand  age  at  time  of  initial  infection  by  dwarf  mistletoe 
may  be  varied  as  desired.  Other  control  variables  include  stand  age 
at  initial  thinning,  stocking  goals,  and  frequency  of  thinning.  Stand 
conditions  and  severity  of  dwarf  mistletoe  infestation  change  with 
time  and  in  response  to  intermediate  cuttings. 

Key  Words:     Stand  yield  tables,  timber  management,  forest 

management,  simulation,  Aroeuthobium  amerioanum, 
Pinus  contovta. 
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Simulating  Yields  of  Managed,  Dwarf 
Mistletoe-Infested  Lodgepole  Pine  Stands 


Clifford  A.  Myers,  Frank  G.  Hawksworth, 
and  James  L.  Stewart 


Introduction 


Procedures  are  available  for  predicting  yields 
in  healthy,  managed  stands  of  lodgepole  pine 
(Pinus  contorta  Dougl.)  in  the  central  Rocky 
Mountains  (Myers  1967).  More  than  half  the 
lodgepole  pine  stands  in  this  area,  however, 
are  infested  with  dwarf  mistletoe,  Arceuthobium 
americanum  Nutt.  ex  Engelm.  (Hawksworth 
1958).  Management  decisions  must  be  made 
for  each  of  these  stands.  Possible  alternatives 
are:  (1)  sanitation  thinning,  (2)  no  treatment 
with  the  hope  that  the  stands  may  'eventually 
become  merchantable,  or  (3)  destruction  and 
regeneration.  Information  on  potential  yields 
is  needed  to  help  forest  managers  evaluate 
these  alternatives.  Procedures  and  a  computer 
program  for  predicting  yields  in  dwarf  mistle- 
toe-infested stands  are  presented  here  to  meet 
this  need.  A  subsequent  publication  will  dis- 
cuss details  of  application  for  pathologists  and 
forest  managers. 

Information  is  available  on  the  effects  of 
A.  americanum  on  growth  rates  of  individual 
lodgepole  pine  trees.  Weir  (1916)  compared 
dominant  heavily  infected  and  healthy  lodge- 
pole pines  in  Washington  and  Idaho.  His  data, 
based  on  50  infected  trees  65  years  old  and  50 
healthy  trees  60  years  old  in  the  same  stands, 
are  summarized  below.  Initially,  trees  that 
became  infected  were  growing  faster  in  diam- 
eter than  those  that  remained  healthy. 


Radial  growth, 

Height  D.b.h.  last  40  years 

(Feet)  (Inches)  (Inches) 

Healthy               48.5  7.8  2.93 

Infected             35.2  6.3  .93 

Infected,  as 
percent  of 

healthy               73  81  32 


Baranyay  and  Safranyik  (1970)  studied  the 
effects  of  dwarf  mistletoe  on  the  growth  of 
lodgepole  pines  in  Alberta,  Canada.  Volume 
differences  were  significant  only  between 
heavily  infected  and  healthy  trees,  and  the 
amount  of  loss  increased  with  age  of  infection. 
Mortality  was  significantly  greater  in  stands 
on  dry  sites  than  on  moist  sites. 

Hawksworth  and  Hinds  (1964)  demonstrated 
the  effects  of  dwarf  mistletoe  in  stands  less 
than  150  years  old.  The  amount  of  damage 
caused  by  the  parasite  was  most  closely  related 
to  time  since  initial  infection,  but  also,  to  a 
lesser  extent,  to  stand  age  when  initial  in- 
fection occurred.  Reduction  in  total  height 
and  d.b.h.  growth  of  dominant  and  codominant 
trees  averaged  about  0.7  percent  per  year  since 
time  of  infection.  In  heavily  infested  stands, 
growth  losses  and  increased  mortality  reduced 
cubic-foot  volume  by  1  to  2  percent  per  year. 
Their  study  dealt  with  past  effects  of  the  para- 
site in  unmanaged  stands;  it  did  not  provide 
data  that  could  be  used  for  direct  yield  pre- 
diction in  managed  stands. 


Data  Collection 

Methods  used  were  similar  to  those  out- 
lined by  Myers  (1966,  1967)  for  healthy  stands, 
but  with  additional  information  on  the  extent 
of  dwarf  mistletoe  infection  in  each  tree. 

Seventy-nine  temporary  plots  were  estab- 
lished in  even-aged  lodgepole  pine  stands  in 
Colorado  and  southern  Wyoming.  Stands  varied 
in  site  index  from  30  to  85,  and  dwarf  mistletoe 
infestation  ranged  from  none  to  very  heavy. 

The  plots  conformed  to  the  following 
standards: 

1.  Were  uniform  in  site  quality,  range  of  tree 
sizes,  and  stand  density  on  and  adjacent 
to  the  plot. 

2.  Varied  in  area  with  average  stand  diameter. 
Most  plots  contained  approximately  150  trees, 
and  ranged  from  0.020  to  0.345  acre. 
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3.  Supported  even-aged,  thinned  or  unthinned 
stands  that  had  not  been  cut  or  otherwise 
disturbed  within  13  years  prior  to  measure- 
ment. The  3  years  beyond  the  10-year  meas- 
urement period  was  to  allow  for  severe  ad- 
justments in  form  that  would  not  be  expected 
in  continually  managed  stands. 

4.  Had  stand  diameters  reasonably  close  to  those 
possible  in  managed  stands  of  equal  age. 
Some  leeway  was  found  to  be  allowable  in 
this  standard,  but  all  extreme  conditions 
were  rejected. 

The  extent  of  infection  in  each  tree  was 
ranked  on  the  42  plots  with  dwarf  mistletoe 
by  a  6-class  mistletoe  rating  system  (Hawks- 
worth  1961).  The  live  crown  was  divided  hori- 
zontally into  thirds. Each  third  was  rated  as: 

0  =  No  mistletoe. 

1  =  Light  —  less  than  half  the  branches 

infected. 

2  =  Heavy  —  more  than  half  the  branches 

infected. 

The  ratings  for  each  third  of  the  crown 
were  added  to  obtain  the  total  for  the  tree. 
The  rating  for  a  tree  can  therefore  vary  from 
0  (no  infection)  to  6  (heavy  infection  in  all 
thirds). 

Measurements  made  on  each  plot  are  in- 
dicated by  entries  on  the  sample  field  forms 
(figs.  1  and  2).  Only  values  needed  to  carry 
out  necessary  computations  are  shown.  Addi- 
tional information  can  be  obtained,  if  desired, 
especially  if  local  experience  indicates  that 
additional  independent  variables  should  be 
tested  in  the  regression  analyses. 

The  plot  description  form  (fig.  1)  provided 
spaces  for  recording  descriptive  material  and 
the  results  of  computations.  Results  that  need 
appear  only  on  computer  printouts  are  given 
to  provide  a  complete  list. 

Field  measurements  of  living  trees  were 
those  shown  in  columns  two  to  seven  of  figure 
2.  Each  tree  was  given  a  temporary  number  by 
stapling  a  numbered  card  to  it.  This  permitted 
most  efficient  use  of  small  crews.  The  record 
of  figure  2  was  completed  one  column  at  a 
time,  yet  all  data  from  any  one  tree  could  still 
be  identified  as  such.  Diameters  of  all  trees 
on  each  plot  were  measured  with  a  diameter 
tape.  Total  heights  of  all  trees  were  measured 
with  either  a  marked  pole  placed  beside  the 
tree,  a  Suunto  hypsometer,  or  an  Abney  level. 
The  record  of  heights  provided  data  for:  (1) 
site  index  determination  (6-8  trees),  (2)  con- 
struction of  a  height-diameter  curve,  and  (3) 
determining  average  height  of  dominant  and 
codominant  trees  (about  20).  Total  ages  were 
determined  from  borings  of  dominant  and  co- 


dominant  trees  at  the  ground  line.  Ages  were 
determined  for  intermediate  and  overtopped 
trees  when  needed  to  confirm  that  the  stand 
was  even-aged,  but  were  not  used  in  com- 
putations. Radial  growth  of  the  wood  at  breast 
height  was  determined  for  each  tree  by  boring 
along  the  best  estimate  of  average  radius. 

Diameters  (outside  bark)  of  trees  that  died 
during  the  previous  10  years  were  recorded. 
Appearance  of  trees  dead  0  to  10  years  was 
determined  by  examination  of  dead  trees  on 
permanent  plots  and  in  stands  for  which  date 
of  thinning  was  known. 


Information  Needed 

Field  data  described  in  the  previous  section 
were  converted  to  volumes  and  other  values 
for  each  plot  and  per  acre.  Basal  area  and 
other  per-acre  values,  average  stand  diameter, 
and  site  index  were  then  used  as  dependent 
and  independent  variables  to  obtain  the  pre- 
diction equations  used  in  program  LPMIST 
(appendix  1).  The  equations,  shown  as 
FORTRAN  statements  in  the  program  listing, 
contain  only  significant  independent  variables. 

Items  computed  and  uses  made  of  them 
in  LPMIST  were  as  follows: 

1.  Present  diameter  outside  bark  was  converted 
to  past  diameter  outside  bark  by  use  of  radial 
wood  growth  and  equations  that  account  for 
bark  growth  (Myers  1964a). 

2.  Past  number  of  trees  equaled  present  number 
of  live  trees  (fig.  2)  plus  any  mortality 
during  the  previous  10  years.  Past  and 
present  numbers  were  raised  to  an  acre  basis. 
Equations  to  estimate  mortality  in  healthy 
(OUT)  and  diseased  (DIE)  stands  were  com- 
puted from  density,  mortality,  and  other  data. 

3.  Present  age  of  the  main  stand  was  the  aver- 
age of  the  ages  of  dominant  and  codominant 
trees  (column  5  of  fig.  2).  Past  age  equaled 
present  age  minus  10  years. 

4.  Average  height  of  dominant  and  codominant 
trees  was  obtained  from  columns  3  and  4 
of  figure  2.  Prediction  equations  for  height 
(HTSO)  were  determined  from  healthy  stands 
with  densities  within  the  range  of  possible 
management  goals.  Data  from  infected  plots 
were  used  to  derive  the  equation  for  re- 
duction in  height  growth  due  to  disease 
(PCT). 

5.  Site  index  was  calculated  by  methods  out- 
lined by  Alexander  (1966).  Because  heavy 
dwarf  mistletoe  infection  reduces  height 
growth,  site  index  on  heavily  infested  plots 
was  estimated  by  measuring  nearby  healthy 
stands  on  comparable  sites. 
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Lodgepole  Pine  Dwarf  Mistletoe  Yield  Study 
Plot  Description  and  Summations 


Lodgepole  Pine  Dwarf  Mistletoe  Yield  Table  Study 
Description  of  Live  Trees 
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6.  Present  basal  area  was  computed  with  the 
diameters  of  column  2  of  figure  2.  Past 
basal  area  was  obtained  from  the  diameters 
of  column  9  of  figure  2  and  diameters  of 
the  dead  trees.  Basal  areas  were  later  used 
in  predictions  of  d.b.h.,  volume,  and  other 
variables. 

7.  Average  stand  diameter  was  the  diameter 
of  the  tree  of  average  basal  area.  This 
definition  of  average  diameter  was  applied 
every  place  this  variable  was  used  or  recorded 
as  a  result.  Initial  average  diameters  and 
other  variables  from  healthy  stands  were 
used  to  obtain  the  equation  for  average 
d.b.h.  after  10  years  (DBHO).  Data  from 
diseased  stands  provided  the  equation  for 
reduction  in  diameter  growth  due  to  dwarf 
mistletoe  (TEM). 

8.  Cubic-foot  and  board-foot  volumes  on  each 
plot  were  computed  with  tree  volume  equa- 
tions of  the  form  V  =  a  +  bD2  H  (Myers 
1964b,  1969).  Plot  volumes  were  raised  to 
volumes  per  acre  before  use  in  other  com- 
putations. Total  cubic  volumes  were  used 
to  compute  equations  for  stand  volume  in 
cubic  feet  (TOTO  and  TOTT).  Total  volumes 
plus  merchantable  cubic-  and  board-foot  vol- 


umes were  used  to  obtain  equations  for  the 
volume  conversion  factors  (FCTR  and  PROD) 
computed  by  subroutine  LPVOL. 
9.  Present  dwarf  mistletoe  rating  of  each  plot 
was  determined  by  averaging  the  individual 
ratings  for  all  live  trees  on  each  plot.  Esti- 
mates of  periodic  increases  in  dwarf  mistle- 
toe rating  (DMR)  of  the  plot  trees  were  ob- 
tained from  results  of  a  previous  study 
(Hawksworth  and  Hinds  1964). 

Several  prediction  equations  are  used  to 
obtain  dwarf  mistletoe  ratings  in  LPMIST. 
One  equation  for  DMR  predicts  the  initial 
rating  if  the  stand  has  never  been  thinned. 
Other  equations  for  DMR  predict  the  current 
rating  as  an  increase  from  a  past  value. 
An  expected  post-thinning  rating  (DMRT) 
is  computed  if  infection  index  is  not  so 
high  (3.0  or  greater)  as  to  make  thinning 
impractical  and  if  a  thinning  from  above 
has  not  already  been  made.  An  equation 
in  subroutine  LPCUT2  then  predicts  the 
percentage  of  trees  to  be  removed  by  thin- 
ning from  above  (REDT)  to  obtain  the  ex- 
pected rating.  For  subsequent  thinnings, 
DMRT  is  computed  from  intensity  of  thin- 
ning and  rating  prior  to  thinning.  Data 
used  in  regression  analysis  to  obtain  these 
equations  came  from  growth  plot  data  and 
actual  and  simulated  thinnings. 
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A  rule  of  thumb  is  available  for  estimating 
DMR  so  that  not  all  trees  need  to  be  rated. 
For  lodgepole  pine  and  A.  americanum: 

log  DMR  =  -0.8814  +  0.0192  P 

where  P  is  the  percentage  of  infected  trees 
in  the  stand.  This  equation  may  be  used  to 
obtain  a  starting  value  for  DMR  when  future 
conditions  of  actual  stands  are  to  be  simulated. 
10.  Thinning  intensity  in  healthy  stands  or 
where  thinning  from  below  is  to  be  done, 
is  based  on  a  relationship  between  d.b.h. 
and  basal  area: 


Average  stand  d.b.h.  Basal  area 

after  thinning  per  acre 
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These  values,  SQFT  in  subroutine 
LPCUT1  and  LPCUT3,  represent  one  possible 
series  of  densities  that  could  be  used  to 
guide  successive  thinnings.  The  growing 
stock  level  shown  above  is  80;  reserve  basal 
area  remains  constant  at  80  square  feet  after 
stand  d.b.h.  reaches  10.0  inches.  Other  stock- 
ing levels  are  named  the  same  way.  For 
example,  level  100  means  that  reserve  basal 
area  will  be  100  square  feet  when  d.b.h. 
is  10.0  inches  or  larger.  Basal  areas  for 
level  100  and  diameters  smaller  than  10.0 
inches  are  obtained  by  multiplying  each  basal 
area  of  level  80  by  the  amount  100/80. 
Values  for  any  stocking  level,  THIN  orDLEV 
in  LPMIST,  are  computed  similarly. 

Equations  for  DBHP  in  subroutine 
LPCUT1  and  LPCUT3  also  describe  the  tabu- 


lated values.  In  this  case,  diameter  is  esti- 
mated when  basal  area  and  the  desired  stock- 
ing level  are  known.  Variables  BREAK  and 
BUST  indicate  points  where  the  relationship 
of  diameter  to  basal  area  has  been  broken 
into  segments  for  convenience  in  regression 
analysis. 

Growing  stock  levels  to  be  left  after 
thinning  from  below  are  indicated  by  assign- 
ing values  to  THIN  and  DLEV  on  data  card 
type  3.  Each  assigned  value  is  a  growing 
stock  level  or  the  basal  area  left  when  d.b.h. 
after  thinning  is  10.0  inches  or  greater. 

11.  Equations  for  DBHE  (used  as  DBHT)  in 
LPCUT1  and  LPCUT3  and  for  ADDHT  in 
the  main  routine  were  derived  from  data 
obtained  in  a  variety  of  thinned  stands. 
Thinning  may  be  simulated  on  a  computer. 
The  sequence  of  operations  is  as  follows: 

a.  Convert  a  series  of  stand  tables  of  actual 
stands  to  1,000  trees  each.  Assign  the 
appropriate  total  height  to  each  tree. 

b.  Compute  average  d.b.h.  and  average 
height  of  each  stand. 

c.  Create  a  set  of  1,000  randomly  arranged 
diameters,  each  with  its  corresponding 
height. 

d.  Create  groups  of  trees,  based  on  per- 
centages of  trees  to  be  retained  (PRET), 
and  tally  the  largest  (thinning from  below) 
or  smallest  (thinning  from  above)  d.b.h. 
in  each  group.  For  example,  if  25  percent 
of  the  trees  are  to  be  retained,  divide  the 
1,000  randomly  arranged  diameters  into 
250  groups  of  four  trees  each.  Tally  the 
largest  or  smallest  diameter  in  each  group 
of  four  trees.  Also  tally  the  corresponding 
height. 

e.  Repeat  step  d  for  various  percentages  of 
retention  and  for  each  stand. 

f.  Compute  the  average  d.b.h.  (DBHE)  and 
average  height  that  results  from  each  per- 
centage retention  (PRET)  in  each  stand. 

g.  Determine  the  difference  between  each 
average  height  and  the  height  of  the  stand 
before  thinning. 

h.  Use  regression  analysis  to  obtain  equa- 
tions that  predict  diameter  after  thinning 
(DBHE)  and  change  in  height  (ADDHT), 
from  diameter  before  thinning  (DBHO)  and 
percentage  of  trees  retained  (PRET). 

i.  Compare  predicted  and  actual  values  of 
DBHE  and  ADDHT,  using  data  from  actual 
stands,  to  insure  that  adequate  predictions 
can  be  made. 

12.  Values  for  AGEO,  DBHO,  and  DENO  on 
data  card  type  3  are  obtained  by  examining 
a  variety  of  young  stands.  For  each  site 
class,  average  d.b.h.  at  various  ages  is  deter- 
mined for  each  of  several  levels  of  stand 
density. 
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Description  of  Program  LPMIST 

Program  LPMIST  consists  of  a  main  pro- 
gram and  four  subroutine  subprograms.  The 
main  program  performs  most  computations  and 
writes  the  yield  tables.  Three  subroutines 
compute  average  stand  d.b.h.  and  stand  density 
after  thinning.  The  fourth  subroutine  com- 
putes factors  that  are  used  in  the  main  program 
to  convert  total  cubic  feet  to  other  units. 

Operations  performed  by  each  routine  are 
identified  by  the  comment  statements  of  the 
source  program  (appendix  1).  Initial  stand 
conditions  and  values  of  several  control  vari- 
ables are  read  into  computer  memory  in  the 
order  and  format  given  in  the  tabulation  of 
order  and  contents  of  the  data  deck.  The  num- 
ber of  yield  tables  computed  and  printed  is 
determined  by  the  values  assigned  NTSTS  on 
card  type  1  and  MIX  on  card  type  2.  NTSTS 
is  the  number  of  sets  of  tables  to  be  produced. 
MIX  is  the  number  of  tables  in  a  set  or  the 
number  of  growing  stock  levels  (DSTY)  tested. 

Operations  are  performed  in  the  following 
sequence: 

1.  Compute  average  height,  basal  area,  volume, 
and  mistletoe  rating  just  prior  to  initial 
thinning. 

2.  Compute  the  current  mistletoe  rating  after 
thinning,  and  execute  the  thinning  if  current 
rating  is  below  3.0.  If  thinning  is  possible, 
subroutines  compute  the  new  stand  density 
and  average  d.b.h.  as  explained  below.  The 
main  program  then  computes  the  new  average 
stand  height. 

3.  Compute  post-thinning  volumes. 

4.  Compute  amounts  removed  by  thinning  and 
data  describing  conditions  before  and  after 
thinning. 

5.  Advance  d.b.h.,  height,  and  mistletoe  rating 
one  prediction  period,  and  compute  new 
volumes.  Mistletoe  rating  is  computed  as  an 
increase  from  a  previous  value  or  as  a  pro- 
jection from  initial  infection,  depending  upon 
whether  or  not  the  stand  has  been  thinned 
since  infection. 

6.  Print  values  appropriate  to  the  stand  age 
if  thinning  is  not  scheduled. 

7.  If  thinning  is  scheduled,  re-thin  by  return 
of  program  control  to  the  operations 
described  in  item  2. 

8.  Repeat  operations  described  in  items  2  to  7, 
inclusive,  until  stand  age  reaches  the  limit 
set  by  ROTA  on  data  card  type  3.  Only 
one  thinning  in  diseased  stands  will  be  from 
above,  and  simulated  by  LPCUT2.  Subse- 
quent thinnings  will  increase  average  d.b.h. 
and  height,  but  by  lesser  amounts  than 


where  the  smaller  trees  make  up  a  very 
large  percentage  of  those  removed.  This 
effect  has  been  observed  in  subsequent  thin- 
nings of  actual  stands,  and  is  simulated 
by  LPCUT3. 
Zero  punches  in  any  data  card  will  cause 
control  to  move  to  the  end  of  the  program, 
a  diagnostic  message  to  be  printed,  and  ter- 
mination of  the  job. 

Clearcutting  at  the  oldest  age  of  interest 
(ROTA)  is  assumed,  because  of  the  serotinous 
nature  of  the  cones  of  most  lodgepole  pines 
in  Colorado  and  Wyoming.  The  program  can 
be  modified  to  show  volume  reserved  for  shelter- 
wood  or  for  seed  trees,  if  desired. 

Subroutines  LPCUT1  and  LPCUT3  compute 
average  stand  d.b.h.  after  thinnings  that  remove 
many  of  the  smaller  trees  and  thus  produce 
increases  in  average  stand  diameter  and  height. 
Successive  percentages  of  trees  to  be  retained 
(PRET)  are  tested  until  the  relationship  be- 
tween d.b.h.,  basal  area,  and  number  of  trees 
is  mathematically  correct  and  d.b.h.  and  basal 
area  agree  with  the  growing  stock  level  specified 
by  THIN  or  DLEV.  Two  major  loops  are  pro- 
vided in  the  subroutines  because  two  equations 
are  needed  for  estimating  post-thinning  d.b.h. 
(DBHE). 

Subroutine  LPCUT2  uses  thinning  standards 
based  on  the  goals  of  sanitation  thinning,  not 
on  THIN  or  DLEV.  The  reduced  infection 
rating  to  be  attained  (DMRT)  is  computed  by 
the  main  program  as  a  function  of  average 
stand  d.b.h.: 


D.b.h.  Rating 
(Inches) 

2  0 

4  0.5 

6  1.0 

8  1.5 

10  2.0 


LPCUT2  then  computes  the  reduction  in  stand 
density  (REDT)  needed  to  attain  this  goal, 
using  d.b.h.  and  rating  just  prior  to  thinning. 
D.b.h.  after  thinning  (DBHT)  can  then  be 
determined  directly  from  the  same  equations  as 
those  for  DBHE  inLPCUTl.  Successive  approxi- 
mations are  unnecessary  because  percentage  of 
trees  to  be  retained  (PRET)  is  known  before 
DBHE  (as  DBHT)  is  computed. 
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Order  and  Contents  of  the  Data  Deck 


Card 
ty_£e_ 


Frequency- 
read 


Variable 
name 


Columns 


Format 


Description  of  variable 


Once 


Each 
test 


NTSTS 


JCYCL 


1-4 


1-4 


14 


14 


Number  of  sets  of  yield  tables  to 
be  produced.     Each  set  is  based 
on  values  of  the  variables  on 
appropriate  data  cards  of  types  2 
and  3 . 

Interval  between  intermediate 
cuts.     A  multiple  of  RINT. 


MIX  5-8  14  Number  of  stocking  levels  or 

values  of  DLEV  to  be  examined  in 
one  test. 


Each 
test 


AGEO 


1-5 


F8.3  Initial  age  to  be  shown  in  a 

yield  table.  Stand  age  when 
first  thinning  occurs. 


DBHO 


9-16  F8.3 


Average  stand  d.b.h.  just  prior 
to  initial  thinning  at  age  AGEO 
and  with  density  DENO. 


DENO 


17-24  F8.3 


Number  of  trees  per  acre  just 
prior  to  initial  thinning  at  age 
AGEO. 


DSTY  25-32  F8.3  Lowest  growing  stock  level  in  a 

test  for  intermediate  cuts  after 
initial  thinning.     Level  will 
increase  by  10  as  many  times  as 
specified  by  MIX. 


RINT  33-40  F8.3 

ROTA  41-48  F8.3 

SITE  49-56  F8.3 

THIN  57-64  F8.3 

START  65-72  F8.3 


Number  of  years  for  which  growth 
and  infection  equations  make  one 
projection  of  growth  or  change. 
Value  is  10.0  for  the  equations 
given  in  Appendix  1. 

Final  age  for  which  stand  data 
are  to  be  given  in  a  yield  table. 

Site  index  on  which  the  set  of 
yield  tables  is  to  be  based. 

Growing  stock  level  for  initial 
thinning  at  age  AGEO.    May  equal 
DLEV . 

Stand  age  at  which  dwarf  mistle- 
toe infection  begins.     Never  enter 
zero.     Enter  number  larger  than 
ROTA  if  infection  will  not  occur 
during  the  rotation. 
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Modifications  for  Other  Species 

Replacement  of  several  statements  will 
modify  the  program  for  other  species  or  regions. 
Replacements  needed  are: 

1.  Statements  for  SQFT,  DBHP,  BREAK,  BUST, 
and  related  computations  that  contain  the 
ratio  of  DLEV  or  THIN  to  80.0,  if  desired. 
This  change  is  needed  if  standards  for  re- 
serve stand  in  LPCUT1  and  LPCUT3  will  be 
different  from  those  shown  in  the  tabulation 
of  the  previous  section. 

2.  Statements  for  TOTO  and  TOTT,  to  make 
cubic  volumes  per  acre  correct  for  the  species 
and  tree  volume  equations  selected. 

3.  Statements  for  FCTR  and  PROD  in  subroutine 
LPVOL  that  are  correct  for  the  species  and 
tree  volume  equations  selected. 

4.  Statements  for  HTSO,  ADDHT,  and  PCT  so 
that  height  growth,  changes  in  height  due 
to  thinning,  and  reductions  in  growth  caused 
by  dwarf  mistletoe  will  be  appropriate  for  the 
species. 

5.  Statements  for  or  that  include  DMR,  DMRT, 
and  REDT;  to  show  correct  relationships 
for  the  host-parasite  interaction  being  simu- 
lated. 

6.  Statement  for  DBHO,  based  on  a  growth 
study  in  healthy  stands  of  the  species  of 
interest  and  a  statement  for  TEM  to  com- 
pute the  effect  of  mistletoe  on  diameter 
growth. 

7.  Statements  for  DBHE  in  subroutine  LPCUT1 
and  LPCUT3  and  for  DBHT  in  LPCUT2 
that  apply  to  the  species  of  interest. 

8.  Statements  that  describe  periodic  losses  in 
numbers  of  trees  in  both  healthy  (OUT) 
and  diseased  (DIE)  stands. 

9.  Table  headings. 

A  Sample  Problem 

The  sample  problem  described  here  pro- 
vides additional  description  of  the  data  deck 
and  of  the  output  (appendix  2).  It  can  also 
serve  as  a  test  problem  to  check  accuracy  of 
punching  of  the  source  deck  and  to  test  com- 
patibility with  local  equipment. 

Assume  a  forest  composed  of  even-aged 
stands  of  lodgepole  pine  that  differ  in  degree 
of  infection  by  dwarf  mistletoe.  Site  indexes 
range  from  30  to  85  (base  100  years)  and  in- 
fection ratings  range  from  0  to  6. 

Problems  to  be  solved  by  the  manager 
of  such  a  forest  include: 

1.  What  growth  can  be  expected  in  healthy 
stands  of  known  site  quality  for  various 


combinations  of  thinning  frequency  and  in- 
tensity? 

2.  How  is  this  growth  affected  by  various  de- 
grees of  dwarf  mistletoe  infection  and  time 
of  initial  sanitation  thinning? 

3.  On  the  basis  of  potential  yields  of  each 
stand,  is  thinning,  replacement,  or  no  treat- 
ment appropriate  for  the  stand  at  this  time? 

4.  Does  each  treatment  decision  appear  appro- 
priate when  impacts  on  other  forest  re- 
sources are  considered? 

This  series  of  questions  cannot  call  forth 
answers  that  contribute  to  good  land  manage- 
ment unless  all  numerical  results  can  be  esti- 
mated to  a  useful  degree  of  accuracy.  Program 
LPMIST  provides  such  answers  for  the  tree 
component  of  the  forest. 

Our  hypothetical  manager  may  decide  to 
compare  yields  of  healthy  stands  with  those 
initially  infected  at  age  10.  Other  variables 
would  remain  constant  for  both  tests  except 
for  stand  conditions  at  initial  thinning  and 
intensity  of  thinning.  His  data  deck  could 
contain  the  following  values: 

NTSTS  -  3,  for  healthy  stands  (test  1),  dis- 
eased and  thinned  at  age  30  (test 
2),  and  diseased  and  first  thinned 
at  age  50  (test  3). 

JCYCL  -  20  years. 

MIX  -  3,  or  3  intensities  of  thinning  will 
be  examined  in  each  test.  For 
brevity,  not  all  tables  will  be  re- 
produced in  appendix  2. 

AGEO  -  30.0  years  for  two  tests  and  50.0 
years  for  the  third. 

DBHO  -  4.5  inches  for  two  tests,  3.0  inches 
for  the  third. 

DENO  -  1000.0  trees  for  two  tests,  2500.0 
trees  for  the  third. 

DSTY  -  80.0,  lowest  subsequent  thinning 
level  of  the  3  to  be  examined. 

RINT    -  10.0  years. 

ROTA   -  130.0  years. 

SITE     -  70.0  feet,  base  100  years. 

THIN    -  120.0,  level  for  initial  thinning. 

START  -  200.0,  10.0,  and  10.0  years  on  type 
3  data  cards  of  test  1,  test  2,  and 
test  3,  respectively.  Any  number 
larger  than  the  value  of  ROTA 
could  replace  the  200.0  shown. 

These  values  will  provide  data  for  com- 
paring differences  in  yields  between  healthy 
and  diseased  stands  and  between  different  types 
of  diseased  stands.  Values  must  be  read  from 
data  cards  assembled  in  the  order:  (1)  type  1, 
(2)  type  2  of  test  1,  (3)  type  3  of  test  1,  (4) 
type  2  of  test  2,  (5)  type  3  of  test  2,  (6)  type 
2  of  test  3,  and  (7)  type  3  of  test  3.  Additional 
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tests  could  be  made  to  examine  the  effect  of 
variations  in  thinning  intensity  or  in  any  other 
control  variable. 

Tables  produced  by  LPMIST  can  be  used 
in  many  ways  to  assist  in  decisionmaking. 
Yields,  number  of  noncommercial  cuts,  number 
of  scheduled  cuts  that  cannot  be  made,  and 
size  of  the  average  tree  are  some  of  the  values 
produced.  Money  yields  and  rates  earned  can 
be  computed  if  necessary  data  on  costs  and 
stumpage  values  are  available.  Stand  ages  at 
culmination  of  mean  annual  increment,  and 
rates  earned  can  help  the  manager  determine 
suitable  rotations  for  his  working  groups. 

A  manager  examining  the  tables  in  appen- 
dix 2,  for  example,  might  reach  the  following 
conclusions: 

1.  A  stand  initially  infested  at  age  10  and  then 
left  without  treatment  for  40  years,  can  pro- 
duce only  a  few  posts  and  some  fuelwood 
by  age  130  (fifth  table  of  appendix  2). 
Planned  frequency  and  intensity  of  thinning 
will  have  no  influence  because  the  stand  is 
already  too  heavily  infested  by  age  50  for 
treatment  to  produce  improvement. 

2.  A  stand  infested  by  dwarf  mistletoe  at  age 
10  but  thinned  at  20-year  intervals  beginning 
at  age  30,  can  produce  useful  wood  products. 
Yields,  including  thinnings,  would  be  less 
than  those  from  healthy  stands  with  the 
same  site  index  and  thinned  according  to 
the  same  schedule.  Also,  actual  yields  of 
diseased  stands  would  be  less  than  the  com- 
puted volumes  because  no  reduction  has 
been  made  for  amounts  of  wood  lost  due  to 
pitch  or  distorted  grain. 

3.  In  healthy  stands,  largest  yields  would  be 
produced  with  relatively  light  thinnings,  such 
as  to  level  100.  In  thinned,  diseased  stands, 
largest  yields  would  be  produced  with  heavier 
thinnings,  such  as  to  level  80.  Comparing 
the  best  yields  in  thinned  stands  with  and 
without  dwarf  mistletoe,  diseased  stands  pro- 
duce about  70  percent  of  the  merchantable 
cubic-foot  and  board-foot  volumes  of  healthy 
stands. 
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APPENDIX  1 


Listing  of  Program  LPMIST 


PROGRAM  LPMIST 
1( INPUT, OUT PUT, TAPE5= INPUT, TAPE6=QUTPUT) 

TO  COMPUTE   AND  PRINT  YIELD  TABLES  FOR   EVEN-AGED   STANDS   OF   LODGE  POLE 
PINE    INFECTED  BY   DWARF  MISTLETOE. 

DEFINITIONS  OF  VARIABLES. 

ADDHT   *    INCREASE   OR   DECREASE    IN  AVERAGE   STAND  HEIGHT   RY  THINNING. 

AGEO  =    INITIAL   AGE    IN   YIELD  TABLE. 

BASC   =   BASAL   AREA  CUT   PER  ACRE. 

BASO  =   BASAL   AREA  PER   ACRE   BEFORE  THINNING. 

BAST  =   BASAL   AREA  PER   ACRE   AFTER  THINNING. 

BDFC  =   BOARD  FEET  CUT   PER  ACRE. 

BDFO  =   BOARO   FEET   PER   ACRE   BEFORE  THINNING. 

BDFT   =   BOARD   FEET   PER   ACRE   AFTER  THINNING. 

CFMC   =   MERCHANTABLE   CU.   FT.   CUT   PER  ACRE. 

CFMO  =  MERCH.   CU.   FT.   PER   ACRE   BEFORE  THINNING. 

CFMT   =  MERCH.   CU.   FT.    PER   ACRE   AFTER  THINNING. 

OBHO   =   AVERAGE   STAND  D.B.H.   BEFORE  THINNING. 

DBHT   =   AVERAGE   STAND  D.B.H.    AFTER  THINNING. 

DENC   =   TREES  CUT   PER  ACRE. 

OENO  =   TREES   PER   ACRE   BEFORE  THINNING. 

DENT   =   TREES   PER   ACRE   AFTER  THINNING. 

DIE   =   TREES   LOST    IN  DISEASED   STANDS   IN    10   YEARS,    IN  PERCENT. 
DLEV   =   GROWING  STOCK   LEVEL   FOR    INTERMEDIATE   CUTS   AFTER   THE  FIRST. 
OMR   =   DWARF  MISTLETOE    INFECTION  RATING. 

DMRT   =  MAXIMUM   INFECTION   EXPECTED   IN   STANDS  AFTER   THINNING.  GOAL 

FOR   STANDS  NOT   ALREADY   BEYOND  OMR  OF  3.0. 
OSTY   =  LOWEST   VALUE   OF  DLEV  USED    IN  A  TEST. 
HTSO  =   TREE  HEIGHT   BEFORE  THINNING. 
HTST   =   TREE  HEIGHT   AFTER  THINNING. 
JCYCL   =    INTERVAL   BETWEEN    INTERMEDIATE  CUTS. 

KSTEP  =    INDICATOR   WITH  VALUE  OF   ONE    IF   CURRENT   THINNING    IS  FROM 
BELOW  AND  TWO   IF  CURRENT   THINNING   IS  FROM   ABOVE . 

KTR  =    INDICATOR   WITH   VALUE   GREATER   THAN   ZERO   IF    A  SCHEOULEO 

THINNING  HAS  BEEN  SKIPPED  BECAUSE  MISTLETOE  INDEX  IS  TOO  HIGI- 
OR   BECAUSE   STAND   IS   ALREADY   TO   SPECIFIED  STOCKING. 

MIX   =  NUMBER   OF    STOCKING  LEVELS  EXAMINED   PER  TEST. 

NFL  AG  =    INDICATOR   WITH  VALUE  GREATER   THAN   ZERO   IF   A   THINNING  FROM 

ABOVE   HAS   BEEN  MADE   AT   ANY  TIME. 
NTSTS   =   NUMBER  OF   TESTS  PER  BATCH. 
OUT   =   PERCENT   MORTALITY   IN  HEALTHY  STANDS. 

PCT   =   PERIODIC   HEIGHT    INCREASE    IN   INFESTEO  STAND,    AS   A  PERCENTAGE 
OF   THE    INCREASE    IN  COMPARABLE   HEALTHY  STANDS. 

PRE T   =   PERCENTAGE   OF   TREES  RETAINED  AFTER  THINNING. 

REDT   =   PERCENTAGE   REDUCTION    IN  NUMBER   OF   TREES   WHEN  DMR  IS 
REDUCED   TO   DMRT   BY  THINNING. 

R INT   =  NUMBER  OF   YEARS  FOR   WHICH  A   SINGLE   PROJECTION   IS  MADE. 

ROTA  =   FINAL   AGE    IN  YIELD   TABLE . 

SITE  =   SITE  INDEX. 

START   =   STANO  AGE  AT   TIME   OF   INITIAL  INFECTION. 

TEM   =   PERIODIC   D.B.H.    INCREASE    IN    INFESTEO   STAND ,    AS   A  PERCENTAGE 

OF   THE    INCREASE    IN  COMPARABLE   HEALTHY  STANDS. 
THIN   =   GROWING   STOCK   LEVEL   FOR    INITIAL  THINNING. 
TOTC   =   TOTAL   CUBIC   FEET  CUT   PER  ACRE. 
TOTO  =  TOTAL   CUBIC   FEET  PER   ACRE   BEFORE  THINNING. 
TOTT   =   TOTAL   CUBIC   FEET  PER   ACRE   AFTER  THINNING. 

COMMON   BA, BAST, DB HO, DBHT, DENO, DMR, DMRT, FCTR,PRET, PROD, REST, VDM 
DIMENSION  VAR ( 9 ) , TEMHI 2) 

READ  NUMBER  OF  TESTS   PER   BATCH  FROM  CARD   TYPE   ONE . 

READ   (5,5)  NTSTS 
5  FORMAT    ( 14) 

IFINTSTS    -LE.   0)   GO   TO  310 

EXECUTE   PROGRAM   ONCE   FOR   EACH   SET  OF    INITIAL   VALUES   OF  INTEREST. 

DO   300    1=1, NTSTS 

READ   INITIAL   VALUES,   ONE   TEST  AT  A   TIME,    FROM  CARD  TYPES   2   AND  3. 

READ   I  5,  10)   JCYCL, MIX 
10   FORMAT  (214) 

IFIJCYCL    .LE.   0   .OR.   MIX   .LE.    0)    GO   TO  310 

READ   (5,15)   AGEO, DBHO, OENO, DSTY.RINT, ROTA, SITE, THIN, START 
15   FORMAT  I9F8.3) 


VARI  1  )  ■■ 

VAR ( 2 )  ■ 

VAR ( 3 )  ■ 

VAR ( 4 )  = 

VARI5I  = 

VAR (6)  s 

VARI7)  = 

VAR ( 8 )  = 

VARI9I  = 
DO  20  L= 


AGEO 
OBHO 
DENO 
DSTY 
R  INT 
ROTA 
SITE 
THIN 
START 
1,9 


IFIVAR(L)  .LE. 
20  CONTINUE 
DLEV   =  0.0 


0.0)   GO   TO  310 


PROVIDE   FOR   SEVERAL   GROWING   STOCK  LEVELS   PER  TEST. 

DO  300  M=l,MIX 
A  =  M 

ADDHT   =  0.0 
BOFO  =  0.0 
BDFT   =  0.0 
CFMO  =  0.0 
CFMT   =  0.0 
DMR  =  0.0 
DMRT   =  0.0 
HTCUM   =  0.0 
KSTEP  =  1 


KTR   =  0 
NFL  AG   =  0 
TIME   =  0.0 

DLEV   =    (DSTY   +    (A   *    10.0))   -  10.0 
BASO   =   DENO  *   0.0054542  *   OBHO  »  DBHO 

COMPUTE  CURRENT   DWARF   MISTLETOE   RATING,   UNTHINNED  STANDS. 

TIME   =   AGEO  -  START 
IFITIME    .LE.   0.01    GO   TO  25 

DMR   =   0.31572   *   0.08654  *   TIME  -   0.00016   *  DENO 
IFIDMR    .LT.   0.0)   DMR   =  0.0 
IFIOMR   .GT.   6.0)   DMR    =  6.0 

OBTAIN  AVERAGE   HEIGHT  AND  VOLUMES  PER  ACRE. 

25    I F ( AGEO   .GT.    45.0)    GO   TO  30 

HTSO  =  3.86111  -  0.05979  *  AGEO  +  0.01215  *  AGEO  *  SITE 
GO   TO  35 

30  HTSO  =  0.33401  -   33.2866   /   AGEO   «•   0.92341   *  ALOGIO(SITE) 
1*   ALOGIO(SITE)    /  AGEO 

HTSO   =    10.0  **  HTSO 
35   PCT   =    1.0  -   0.0165   «   DMR  *  DMR 

HTSO   =   HTSO   *  PCT 

COMPUTE   TOTAL   CU.    FT.   ANO  CONVERT   TO  OTHER  UNITS. 

D2H  =   DBHO   *   DBHO  *  HTSO 
IFID2H   .GT.    7000.0)   GO   TO  40 

TOTO  =  (0.00276  *  D2H  -  0.00059  *  BASO  -  0.00577)  *  DENO 
GO  TO  45 

40   TOTO  =    (0.00248   *   D2H   ♦    1.96336)    *  DENO 
45    I F ( DBHO    .LT.    5.0)   GO   TO  50 

VDM   =  DBHO 

BA  =  BASO 

CALL  LPVOL 

BOFO   =   TOTO   *  PROD 

CFMO   =   TOTO   *   FC  TR 
50   REST   =  THIN 

ENTER  LOOP   FOR  REMAINING  COMPUTATIONS  AND  PRINTOUT. 

DO  250  K=l, 100 

I F ( AGEO    .GE.   ROTA )   GO   TO  90 

COMPUTE   D.B.H.   AFTER  THINNING. 

IFIDMR    .LT.    3.0)    GO  TO  55 

BAST  =  BASO 
DBHT  =  DBHO 
DMRT   =  DMR 
HTST   =  HTSO 
KTR  =  1 
GO   TO  75 
55    IFIDMR   .EO.   O.lM    Gn   TO  61 
IFINFLAG    .GY.   01    G(l   TO  62 
OMRT   =   0.25   t   DPHO  -  0.5C 


0.(1)  DMRT  =  0.0 
RMR  )    r,0   TO  62 


0.0279   *  PRHT 


"ASOI    GO   TO  70 


IF (DMRT  .LT. 
IFIOMRT  .GF. 
CALL  LPCUT2 
NFLAG   =  1 
KSTEP   =  2 
GO  TO  65 

62  CALL  LPCUT3 
KSTEP   =  I 
DMRT   =   OMR  * 
GC   TO  65 

63  DURT   =  DYK 
CALL  LPCUT1 
KSTFP  =  1 

65  IFIBAST  .LT. 
BAST  =  BASO 
DBHT  =  0I1HQ 
DMRT  =  DMR 
HTST  =  HTSO 
KTR  =  1 
GO  TO  75 

C 

C   COMPUTE   HPIGHT   AMD   VOLUMES   AFTEx  THINNING. 

C 

70  GO   TO    17  1,72),   KSTC  P 

71  ADDHT   =   6.  79950  -   3.41979   <■   ALOG  I  0  ( PR  E  T ) 
GO   TU  73 

72  AOOHT   =   3.76362   *   ALOG  1 0  (  PRF.  T  )   -  7.97347 

73  HTCUM  =   HTCUM   ♦  ARMHT 
HTST  =   HTSC   *  ADDHT 

75   JDENT    =    I  HAS  T   /    (0.0054542   «   DBHT   *   DBHT))    ♦  0.5 
DENT   =  JDENT 

BAST   =   0.0054542   *   D»FT   *   DhHT   *  HENT 
D2H  =   DBHT   *   DBHT   *  HTST 
IFID2H   .GT.    7C00.0)    GC   TO  80 

TCTT  =  (0.C0276  *  D2H  -  0.00059  s  ("AST  -  C.C0577)  *  DENT 
GO  TO  85 

80   TCTT   =    (0.OO24R   «   02 H   ♦    1.96336)    *  DENT 


CONVERT    TOTAL   CU.    FT.    TO  OTHER  UNITS. 

85   I F ( DBHT    .LT.    5.0)    GO  TO  90 
VDM   =  DBHT 
RA  =  BAST 
CALL  LPVOL 
BOFT   =   TCTT   *  PROD 
CFMT   =    TOTT   *   FC  TR 

CHANGE   MODE   AND   ROUND  OFF   FOR  PRINTING. 
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90  J AGEO 
JSITE 
JOENO 
JHTSO 
JTOTO 
JTOTO 
JBASO 
JCFMO 
JCFMO 
JBOFO 
J  BDFO 
JHTST 
J  TOT  T 
J  TOT  T 
JCFMT 
JCFMT 
IF( JCFMT 
J  BOFT  = 
JBDFT  = 
[F( JBDFT 
JBAST 
JDENC 
JBASC 
JTOTC 
JCFMC 
IFI JCFMC 
JBOFC 
I F ( J  BOFC 


AGEO 
SITE 

DENO   *  0.5 
HTSO   *  0.5 
(  TOTO   »   0. I) 
JTOTO   *  10 
BASO  ♦  0.5 
(CFMO  *  0.11 
JCFMO  *  10 
( BDFO   *  0.01) 
JBDFO   *  100 
HTST   +  0.5 
(  TOT  T 
JTOTT 
(  CFMT 
JCFMT 
.  GT  . 

0.01  ) 
100 
.GT.  JBDFOI 
BAST   +  0.5 
JOENO  -  JDENT 


0.  5 


(  BDFT 
JBOFT 


*  0.  1  ) 

*  10 

*  0.1)    ♦  0.5 

*  10 

JCFMO)    JCFMO   =  JCFMT 
.5 


JBDFO   =  JBOFT 


JBASO 
JTOTO 
JCFMO 

LE. 
JBOFO 

LE. 


JBAST 
JTOTT 
JCFMT 

)  JCFMC 
JBDFT 

I  JBDFC 


WRITE  HEADINGS   FOR   YIELD  TABLE. 

IF(K    . GE •   2)    GO   TO  120 
WRITE   (6,95)    JSITE, THIN, DLEV 
95   FORMAT    (  1 H 1 , / / / , 39X , 5 3H Y I  ELDS  PER   ACRE   OF   EVEN-AGED   STANDS   OF  LODG 
1EP0LE   PINE/1H   ,57X,11HSITE    INDEX    ,I3/1H    , 38 X , 29HTH I NN I NG  INTENSITY 
2-   INITIAL-   ,F5.0,2X, 12HSUBSE0UENT-  ,F5.0) 
WRITE  16,100) 

100   FORMAT    I IHO, 25X, 38HENTIRE   STAND   BEFORE   AND   AFTER    TH I NN I NG , 28X , 26HP 
1ERI0DIC    INTERMEDIATE   CUTS ) 
WRITE  (6,1051 

105   FORMAT    ( 1H0.9X.5HST AN D,10X,5H BASAL, 3X.7H AVE RAGE, 2X.7HAVER AG E ,3X, 5H 
1T0TAL , 3X , 9HMER CHANT-, 3X, 9HSAWT I MBE R , 9X , 5HBASAL,4X,5HT0TAL,3X,9HMER 
2CHANT-,3X,9HSAWTIMBER) 
WRITE  (6,110) 

110   FORMAT    (  1H   ,  I0X , 3HAGE , 4X , 5H TREE S , 3 X , 4HARF A , 4X , 6H0 . B . H . , 3 X , 6HHE I GHT 
1,2X,6HV0LUME,2X,11HABLE   VOL UME, 4X.6H VOLUME,  3X,5HTREES,3X,4HAREA,3X 
2,6HV0LUME,2X, 11HABLE   VOLUME , 4X , 6HV0LUME ) 
WRITE  (6,115) 

115  FORMAT  (  1H  ,  8X ,  7H  (  YEAR  S  )  ,  3X  ,  3HN0.  ,  3  X  ,  6HS0  .  F  T.  ,  4  X  ,  3H I  N.  ,  6  X  ,  3HFT  .  ,  4X 
1 , 6HCU.FT. , 5X, 6HCU.FT. , 6X,6HBD.FT.,4X,3HN0.,3X,6HSQ.FT.,2X,6HCU.FT. 
2,5X,6HCU.FT.,6X,6HBD.FT.) 

WRITE   TABLE   ENTRIES   OF   DIAMETER,    VOLUMES,  ETC. 

120  WRITE    (6,125)    JAGEO,JDENO, JBASO, DBHO, JHTSO, JTOTO, JCFMO, JBDFO 
125   FORMAT  I1H0,9X,I4,4X,I5,2X,I4,5X,F5.1,5X,13,4X,15,6X,I5,6X,I6) 
I F ( AGEO   .GE.   ROTA )   GO   TO  255 

WRITE    I  6,1 30)   JAGEO, JDENT, JBAST, DBHT, JHTST, JTOTT, JCFMT, JBDFT, JDENC 
1 , JBASC, JTOTC, JCFMC, JBDFC 
130   FORMAT    (  1H   ,  9X ,  I  4 ,  4X ,  I  5  ,  2  X  ,  I  4 ,  5  X  ,  F  5  .  1  ,  5  X  ,  I  3  ,  4X  ,  I  5  ,  6  X  ,  I  5  ,  6X  ,  I  6  ,  4X  ,  I 
15,  3X, I  3, 5X, 14, 6X, 14, 8X, I  5 ) 


COMPUTE   VALUES   FOR   EACH  PERIOD.    THIN   AS  SPECIFIED. 

I R 1  NT   =   R I  NT 

IK   =   JCYCL   /  IRINT 

UO  200  L=1,IK 

AGED  =   AGEO   ♦  R I  NT 

I F ( AGEO    .GT.   ROTA)    GO   TO  255 

COMPUTE   CURRENT   DWARF   MISTLETOE  RATING. 


TIME   =   AGEO  -  START 
I F ( OMR    .GT.   0.0)   GO   TO  135 
IFITIME    .LE.    0.0)   GO   TO  150 
OMR   =  0.31572   +   0. 08654   *  TIME 
GO   TO  145 

1.0)    GO   TO  140 
0.07   *  RINT 


0.00016   *  DENT 


(0.03   +   0.038  *   DMRT)    4  RINT 


135    I F ( DMRT  .LE 
OMR  =  DMRT 
GO   TO  145 
140   DMR   =  DMRT 

IFI L    .LE.   2)    GO  TO  145 
DMR  =   DMR   ♦  0.07   ♦  RINT 
145    I F I DMR    .LT.   0.0)   DMR   =  0.0 
I F I DMR    .GT.    6.01    OMR   =  6.0 

COMPUTE   NEW   D.B.H.   BEFORE    THINNING  AND  ROUND   OFF   TO  0.1  INCH. 

150   DBHO  =    1.0222*DBHT   ♦   0.0151*SITE   -    1 . 2 4 1 7*AL0G 1 0  I B AS T )    +  2.1450 
I  Ft  DMRT    .LE.    3.9)   GO   TO  155 

T EM   =    I DBHO  -  DBHT)    *    11.0  -    (0.192   •>   DMRT  -  0.754)1 
DBHO   =   DBHT   +  TEM 
155    I DBHO  =   DBHO   *    10.0   +  0.5 
DBHO   =  IOBHO 
DBHO   =   DBHO  *  0.1 
I F ( DENT    .GT.    1000.0)   GO   TO  160 
DIE   =    (3.81   *  DMRT  -   6.63)    *  0.01 
IFIDIE    .LT.   0.0)    DIE   =  0.0 
GO  TO  165 

160   DIE   =    18.64  ♦   3.28  »   DMRT)    *  0.01 
165  OUT   =  0.0 

IFIOBHT    .GE.    10. 0)   GO   TO  170 

OUT   =   0.05285  -   0.01346  *   DBHT   t   0.00226   »   DBHT   *   OBHT   +  0.0000066 


1»   BAST   *   BAST   -   0.0001931   *   OBHT   *  BAST 
IFIOUT    .LT.   0.0)   OUT   =  0.0 
170   IFIDIE   -LT.   OUT )   OIE   =  OUT 

J  DENO   =    (DENT   *    (1.0  -   DIE))    *  0.5 
DENO   =  JDENO 

BASO  =   OENO  *    (0.0054542   *   DBHO   *   DBHO ) 

OBTAIN  AVERAGE  HEIGHT  AND   VOLUMES  PER  ACRE. 

DO   180  J=l,2 
LUB   =  J 

GO  TO   (  172,  174) ,  LUB 
172  YARS   =  AGEO 

GO  TO  176 
174  YARS  =   AGEO  -  RINT 
176    IFIYARS    .GT.    45.0)   GO   TO  178 

TEMH(J)    =   3.86111   -   0.05979  *   YARS  +   0.01215   *  YARS   *  SITE 

GO  TO  180 

178   TEMHIJ)    =   0.  33401   -    33.2866   /   YARS   <■   0.92341   *  AL0G10ISITE)    +  6.2 
1811   *   AL0G10ISITE)    /  YARS 

TEMHIJ)    =    10.0   *»  TEMHIJ) 
180  CONTINUE 

PCT   =    1.0  -   0.0028  *   DMRT   *   DMRT   *  DMRT 

CHNG  =    (TEMH(l)    -   TEMH(2>)   *  PCT 

HTSO   =   HTST  +  CHNG 

COMPUTE   TOTAL   CU.    FT.   AND  CONVERT   TO  OTHER  UNITS. 

02H   =   DBHO  *   DBHO  *  HTSO 
IFID2H   .GT.    7000.0)   GO   TO  185 
TOTO  =    (0.00276   *   D2H  -   0.00059   *   BASO  - 
GO  TO  190 

185   TOTO  =    (0.00248  *   D2H  -f    1.96336)    *  DENO 


0.00577)    *  DENO 


190    I F ( DBHO  .LT. 
VOM   =  DBHO 
BA  =  BASO 
CALL  LPVOL 
BDFO  =  TOTO 
CFMO  =  TOTO 


5.0)   GO   TO  195 


PROD 
FCTR 


CHANGE   MODE   AND  ROUND  OFF   FOR  PRINTING. 


195    IFIL    .EO.    IK)    GO   TO  205 


KDENO 

DENO   +  0.5 

K  AGEO 

AGEO 

KHTSO 

HTSO   t  0.5 

KBASO 

BASO   +  0.5 

KTOTO 

(TOTO  *   0.11    ♦  0.5 

KTOTO 

KTOTO   *  10 

KCFMO 

( CFMO   *   0.1)    +  0.5 

KCFMO 

KCFMO  *  10 

KBDFO 

(BDFO  *   0.011    ♦  0.5 

K80F0 

KBDFO   »  100 

WRITE  VALUES   FOR   THE   PERIOD   IF    THINNING   IS   NOT  DUE. 


,  F4 


WRITE  (6,125)  KAGEO, KDENO, KBASO, DBHO, KHTSO, KTOTO, KCFMO, KBDFO 
OBHT  =  DBHO 
BAST  =  BASO 
OENT  =  OENO 
DMRT  =  DMR 
HTST  =  HTSO 
200  CONTINUE 

PREPARE   TO   START   LOOP  AGAIN   FOR   NEXT  THINNING. 

205  REST   =  DLEV 
250  CONTINUE 

255    I F I  START    .GE.   ROTA)   GO   TO  265 

WRITE   (6,260)    ST AR T , DMR , ROTA 
260   FORMAT    (  1  HO, 25X , 4 1H0WARF   MISTLETOE    INFECTION   STARTED  AT  AGE 
116H   ANO  RATING  WAS   ,F5.1,8H   AT   AGE  ,F4.0I 

GO  TO  275 
265   WRITE    (6,270)  ROTA 
270   FORMAT    I 1H0, 25X, 63HDWARF   MISTLETOE    INFECTION  DID  NOT   OCCUR  OURIN 

1THE  ROTATION  OF   ,F4.0,7H   YEARS. ) 
275   1FIKTR    . EO.   0)   GO   TO  285 

WRITE  16,280! 

280   FORMAT    I  IHO ,  25X ,  52HN0TE   THAT  NOT   ALL    SCHEDULEO   THINNINGS   WERE  PO 

1  ISLE. ) 
285  WRITE  (6,290) 

290   FORMAT    I  1  HO , 25X ,  66HMERCH.   CU.    FT.    -   TREES   6.0   INCHES   D.B.H.  AND 
1RGER   TO   4-INCH   TOP. I 
WRITE  (6,295) 

295   FORMAT    I IHO, 25X, 59HB0.    FT.   -   TREES   6.5    INCHES  O.B.H.    AND  LARGER 
1   6-INCH   TOP . ) 

PREPARE   FOR  NEXT   TABLE   OF   THE  TEST. 

AGEO   =   VARI 1 ) 
DBHO  =   VARI 2) 
DENO   =   VARI 3 ) 
300  CONTINUE 
GO  TO  350 

PROGRAM  CONTROL   GOES   HERE    IF   ANY   ZEROS   IN   DATA  DECK. 
3L0  WRITE   (6, 3201 

320   FORMAT   ( IH 1 , / // ,  10X , 64HE X  EC UT I  ON   STOPPED   BECAUSE   OF   NEGATIVE  OR 

1R0    ITEM  ON  DATA   CARD.  ) 
350  CALL  EXIT 
END 
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SUBROUTINE  LPVOL 
TO  CONVERT   TOTAL   CU.   FT.    TO  MERCH. 


CU.   FT.    AND   TO  BO.  FT. 


COMMON  BA, BAST, DBHO, DBHT , OENO , OMR, DMR  T  ,  FC  TR  ,  PRE  T  , PROD, RES T.VDM 
FCTR   *  0.0 
PROO  -  0.0 

[FIVOM   .LT.    5.0)   GO  TO  10 


OBTAIN  CONVERSION  FACTORS  FOR  MERCH.  CU.  FT. 
IN  TREES   6.0   INCHES  O.B.H.    AND  LARGER. 


VOLUMES   TO  4.0-INCH  TOP 


IFIVOM   .GT.   6.7)  GO  TO  2 

FCTR   =   0.31963  »  VDM  -  1.42291 
GO  TO  6 

2    IFIVOM   .GT.   9.8)  GO  TO  4 

FCTR  =  3.68255  -  0.14007  »  VOM  -  13.54644  /  VOM 
GO  TO  6 

4  FCTR   =   0.99503  -  0.58018   /  VDM 

6   IFIVOM   .LT.   8.0)  GO   TO  10 


OBTAIN  CONVERSION  FACTORS  FOR  BO.  FT. 
6.5    INCHES   O.B.H.   AND  LARGER. 


VOLUMES   TO  6-INCH   TOP    IN  TREES 


IFIVDM   .GT.    10.0)   GO  TO  8 

PROO  =  2.08874  ♦  0.18091  *  VDM  ♦  0.00045  *  BA 
GO  TO  10 

8   PROD  =  0.16583  ♦   3.74174   *  AL0G10IV0M) 
10  RETURN 
END 

SUBROUTINE  LPCUT1 


ALOG 10  I  PRE  T ) 
ALOG 1 0  I DBHO ) 


♦   0.79921   *   ALOG 1 0  I DBHO ) 


0.00816  *  (PRET 


(PRET  -  50.0) 
♦  0.5 


0.01  ) 


DBHE 


/  80.0 
GO   TO  50 

(0.08682 


19  *    (PRET  -   50.0)  * 
IOBHE   =   DBHE   *  10.0 
DBHE  =   1 DBHE 
DBHE  =  DBHE  *  0.1 
DENE   =   OENO  *  (PRET 
NOENE  -  DENE  +  0.5 
DENE  =  NDENE 
BASE  =   0.0054542   *  DBHE 
NBASE   =   BASE  *   10.0  +  0 
BASE   =  NBASE 
BASE  =  BASE  *  0.1 
BREAK   =   49.9  *  REST 
IFIBASE   .GT.   BREAK ) 
DBHP   =    (80.0   /  REST) 
GO   TO  52 

BUST  =  66.2  *   ( REST   /  80.0) 
IFIBASE    .GT.    BUST)   GO   TO  51 
DBHP  =    180.0   /  REST)   *  (0.10938 
GO  TO  52 

TMPV  =  BASE  *   (80.0   /  REST) 
TEM  =   TMPY  *  TMPV 

DBHP   =    19.04740  *   TMPY  -  0.26673 
1   -  448.76833 
I F ( TMPY   .GT.   80.0)   DBHP   =  DBHO  ♦ 
I DBHP  =   DBHP  *   10.0  ♦  0.5 
DBHP   =    I DBHP 
DBHP   =   DBHP   *  0.1 
IF! DBHP  -   DBHE)  60,70,61 


*   BASE)    +  0.94636 


*   TEM  +   0.0012539   *   TEM  *  TMPY 


TO  ESTIMATE    INCREASE    IN  AVERAGE   O.B.H.   OUE   TO  THINNING   LODGE  POLE  PINE 
IF  DWARF  MISTLETOE  RATING   EOUALS  ZERO. 

COMMON   BA, BAST, 0BHO,0BHT, OENO, DMR, DMRT, FCTR, PRET, PROO, REST, VDM 
IFIOBHO   .LT.    9.5)   GO   TO  30 

COMPUTE   O.B.H.    IF   DBHO    IS   LARGE   ENOUGH  FOR   BASAL   AREA   TO  REMAIN  CONSTANT. £ 

PRET  =    100. 0 
00  21  KJ=1,100 
I F I PRET   .LT.    50.0)   GO   TO  5 

DBHE  =   0.44222  ♦   1.03170  *  DBHO  -   0.00816  *    (PRET   -   50.0)    -  0.0000 
19   *   (PRET  -   50.01    *   (PRET  -  50.0) 
GO  TO  11 

5  PDBHE    -  0.37321  -   0.17274  *   AL0G10 I PRET)    +  0.79921   *  AL0G10 I DBHO) 
1   +  0.09315  *  AL0G10I PRET )    *   ALOG 1 0  I DBHO ) 
DBHE  =   10.0  **  PDBHE 
11    I DBHE  =   DBHE     *    10.0  ♦  0.5 
DBHE  =  IOBHE 
DBHE  =   OBHE  *   0. 1 
DENE   =   DENO  *  PRET   *  0.01 
NDENE  =   DENE   »  0.5 
DENE   =  NOENE 

BASE   =  0.0054542   *   DBHE   *   DBHE   *  DENE 

NBASE  =  BASE  *  10.0  +  0.5  C 
BASE  =  NBASE  C 
BASE  =  BASE  *  0.1  C 
TMPY  =  0.0054542  *  DBHE  *  DBHE  C 
TEM  =  BASE  -  REST 
IFITEM   .LE.   TMPY)   GO   TO  70 

I F I TEM  .LT.  4.0)  GO  TO  20  C 
PRET  =  PRET  -  1.0  C 
GO  TO  21  C 

20  PRET  =  PRET  -  0.3 

21  CONTINUE 
GO  TO  70 

:,    COMPUTE   D.B.H.    IF   BASAL  AREA   INCREASES  WITH  O.B.H. 

30  PRET   =  40.0 

I F I DBHO   .GT.    7.0)   PRET   -  70.0 
DO  65  J=l,  100 

IF! PRET   .GE.   50.0)   GO   TO  40 
PDBHE   =   0.37321  -  0.17274 
1  ♦  0.09315  *  AL0G10I PRET ) 
DBHE   =   10.0  **  PDBHE 
GO  TO  45 

40   DBHE   =   0.44222   +    1.03170  *  DBHO 


60  PRET   =   PRET  *  1.02 
GO  TO  65 

61  PRET   =   PRET   *  0.98 
65  CONTINUE 

70   OBHT   =  DBHE 

COMPUTE   POST-THINNING  BASAL  AREA. 

I F ( DBHT    .GT.    5.0)   GO   TO  75 

SOFT   =    11.58495   »   OBHT  -  11.09724 

GO  TO  76 

75  I F ( OBHT    .GE.    10.0)   GO   TO  77 
TEM  =   DBHT  *  DBHT 

SOFT   =   7.76226  »   DBHT   +0.85289   *   TEM  -0.07952   *   TEM  *  DBHT-3. 45624 

76  BAST   =    (REST   /   80.0)   *  SOFT 
GO  TO  80 

77  BAST   =  REST 
80  RETURN 

END 


SUBROUTINE  LPCUT2 

TO  ESTIMATE  CHANGE   IN  AVERAGE   D.B.H.    DUE   TO   THINNING  LODGEPOLE  PINE 
IF  DWARF   MISTLETOE  RATING  DETERMINES   THE  STANDARDS. 

COMMON   B A, BAST, DBHO, DBHT, DENO, DMR, DMRT, FCTR, PRET, PROD, RE ST, VDM 

COMPUTE   STAND  DENSITY   AFTER   A   THINNING   THAT  REDUCES   THE  INDEX. 

IFIDMR   .LT.    2.0)    GO   TO  5 

R EOT   =   77.5  -   8.5  *   DBHO   +    10.0  *  DMR 

GO  TO  10 

5  REDT   =    15.5  -  8.5  *  DBHO   +   41.0  *  DMR 
10   PRET  =    100.0  -  REDT 

DENT  =   DENO  *   (PRET   *  0.01) 
IOENT   =   DENT   ♦  0.5 
DENT   =  IDENT 

COMPUTE   O.B.H.   AFTER   THINNING  TO  DESIRED  DENSITY. 
I F I PRET    .LT.    50.0)   GO  TO  15 

DBHT   =   0.96559  *   DBHO  +   0.00668   *    (PRET  -   50.0)    ♦   0.00015   *  [PRET 
1   -   50.0)   *   (PRET  -   50.0)   -  0.50568 
GO  TO  20 

15   DBHT   =   0.33478  *   ALOGIO(PRET)   +    1.42477  *   AL0G1 0 ( DBHO )    -  0.21199  * 
1AL0G101 PRET )    *  ALOG10IOBHO)   -  0.67651 

OBHT   =   10.0  **  OBHT 
20   I DBHT   =   DBHT   *   10.0  ♦  0.5 

OBHT   =    I DBHT 

OBHT   =   DBHT   *  0.1 

BAST  =   0.0054542   *   DBHT  *  DBHT  *  DENT 

RETURN 

END 


SUBROUTINE  LPCUT3 

TO  ESTIMATE    INCREASE    IN   AVERAGE   D.B.H.   DUE   TO  THINNING  FROM  BELOW  IF 
DWARF  MISTLETOE  RATING    IS  GREATER   THAN  ZERO. 

COMMON   B A, BAST, DBHO, OBHT, DENO, DMR, DMRT, FCTR, PRET, PROD, REST, VDM 
IFIOBHO   .LT.   9.5)   GO   TO  30 

COMPUTE   O.B.H.    IF  DBHO    IS   LARGE   ENOUGH  FOR   BASAL  AREA   TO  REMAIN  CONSTANT. 

PRET   =  100.0 

DO   21  KJ=1,100 

IFIPRET    .LT.    50.0)   GO   TO  5 

DBHE  =   0.44222   +   1.03170   *  DBHO  -   0.00816   *    (PRET   -  50.0)   -  0.0000 
19  *   (PRET  -   50.0)    *    (PRET  -  50.0) 
GO  TO  11 

5   PDBHE   =   0.37321  -   0.  17274  »   ALOG 1 0 1 PRET )    +  0.79921   *   ALOG 10 1 DBHO ) 
1   *  0.09315  *  AL0G10IPRET)    *  ALOG10IOBHO) 
DBHE  =   10.0  **  PDBHE 
11   TEM  =  DBHE  -  DBHO 

DBHE   =   DBHO   +  TEM   *  0.5 

I  DBHE  =  DBHE  *   10.0  +  0.5 

OBHE  =    I DBHE 

DBHE   =   DBHE   *  0.1 

DENE  =   DENO  *  PRET   *  0.01 

NDENE   =   DENE   +  0.5 

DENE  =  NDENE 

BASE  =   0.0054542   *   DBHE   *  OBHE   *  DENE 

NBASE   =   BASE   «   10.0  +  0.5 

BASE   =  NBASE 

BASE   =   BASE  *   0. 1 

TMPY  =   0.0054542  *   DBHE  *  DBHE 

TEM  =   BASE  -  REST 

IFITEM   .LE.   TMPY)   GO   TO  70 

IFITEM    .LT.   4.0)   GO  TO  20 

PRET   =   PRET  -  1.0 

GO  TO  21 

20  PRET   =■   PRET  -  0.3 

21  CONTINUE 
GO  TO  70 

COMPUTE  D.B.H.    IF   BASAL   AREA   INCREASES  WITH  O.B.H. 

30   PRET   =  40.0 

IFIOBHO   .GT.   7.0)    PRE  T  =  70.0 
DO  65  J=l, 100 

IFIPRET   .GE.   50.0)   GO  TO  40 

POBHE  =  0.37321  -   0.17274  *   AL0G10IPRETI    +   0.79921   *   ALOG 10 ( DBHO ) 
1   +  0.09315  *   AL0G10IPRET)    *   ALOG 1 0 1 DBHO ) 
DBHE  =   10.0  **  PDBHE 
GO  TO  45 

40   DBHE  =   0.44222   »   1.03170  *   DBHO  -   0.00816  *    (PRET  -  50.0)    -  0.0000 
19  *   (PRET  -   50.0)   *    I  PRE  T  -  50.0) 


11- 


45  T EM   =  DBHE 


IFiTMPY    .GT.    80.0)   OBHP   »  DBHO 


DBHE  = 

UdHU    +    ten    *    U  mJ 

52 

i  nnuo    —    i'iijui)          i  n    n    j.    (\  t; 
1 UOMr    —    UHMK    »    1U*U    +    U  •  D 

I OBHE 

—   DBHE  *   10«0  ♦■  0.5 

noun     —      i  nnuD 

Uonr    —  lUtJnr 

DBHE  = 

I  DB  HE 

DBHP   =   DBHP   *   0» 1 

DBHE  = 

OBHE   *  Om 1 

I F ( DBHP   -   DBHE)  60,70,61 

rtChin    *    i  ddct    *    n  nit 

60 

PRET   =   PRET   *  1.02 

NDENE 

=  DENE  +  0.5 

GO  TO  65 

DENE  = 

NDENE 

61 

PRET   =   PRET   *  0.98 

BASE  = 

0.005*542  »  OBHE  *  DBHE  * 

65 

CONTINUE 

NBASE 

=   BASE   *    10.0  +  0.5 

70 

DBHT   =  DBHE 

BASE  - 

NBASE 

C 

BASE  = 

BASE  *   0.  1 

COMPUTE   POST-THINNING  BASAL  AREA. 

BREAK 

=  49.9  *  REST   /  80.0 

I F  t  BASE    .GT.   BREAK )   GO   TO  50 

IF1DBHT   .GT.   5.0)   GO  TO  75 

DBHP  = 

(  80.0   /  REST)    *  (0.08682 

SQFT  =    11.58495  *  DBHT  -  11.09724 

GO  TO 

52 

GO  TO  76 

BUST  = 

66.2  *    (REST   /  80.0) 

75 

I F ( DBHT    . GE •    10.0)    GO   TO  77 

I F ( BASE   .GT.    BUST )    GO   TO  51 

TEM  =   DBHT   *  OBHT 

DBHP  = 

(80.0   /  REST)   *  (0.10938 

*  BASE)    -  0.17858 

SQFT   =   7.76226  *   DBHT   *0. 85289  * 

GO  TO 

52 

76 

BAST   =    (REST   /   80.0)    *  SOFT 

TMPY  = 

BASE  *   (80.0   /  REST) 

GO  TO  80 

TEM  = 

TMPY  *  TMPY 

77 

BAST   =  REST 

DBHP  = 

19.04740  *   TMPY  -  0.26673 

*   TEM  *   0.0012539   *  TFM   *  TMPY 

80 

RETURN 

-  448 

.76833 

END 

APPENDIX  2 
Output  of  Sample  Problem 


YIELDS   PER   ACRE   OF   EVEN-AGED   STANDS  OF   LODGEPOLE  PINE 

SITE    INDEX  70 
THINNING    INTENSITY-    INITIAL-      120.      SUBSEQUENT-  80. 


S  IAND 
AGE 
(YEARS) 

30 
30 


50 
60 


70 
70 


90 
90 


IIP 
110 


120 
130 


TREES 
NC. 


moo 

505 


500 
218 


214 
136 


136 
136 


84 
56 


56 
56 


ENTIRE    STAND   BEFORE   AND   AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

110 

72 


130 
70 


112 

80 


98 
116 


109 

80 


92 
105 


AVEKAGF 
D.8.H. 
IN. 

4.5 
5.1 


6.9 
7.7 


AVERAGE 
HEIGHT 


28 
29 


9.8 
10.4 


12.5 
13.2 


IS. 4 
16.2 


17.4 
18.5 


36 
41 


56 
56 


70 
73 


73 


TOTAL 
VOLUME 
CU.FT. 

1470 
1010 


2660 
15C0 


3140 
2280 


3710 
2560 


3760 
2780 


3300 
3830 


MERCHANT- 
ABLE VOLUME 
CU.FT. 

210 
210 


2000 
1270 


2910 
2140 


3520 
2440 


3600 
2670 


3180 
3690 


SAW  T I M8ER 
VOLUME 
BD.FT. 

0 
0 

0 

0 


12300 
9C0O 


15900 
1  1200 


17300 
13000 


15900 
18800 


TREES 
NO. 


PERIODIC    INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


DWARF   MISTLETCE    INFECTION   Din   MOT   OCCUR   DURING   THE   ROTATION  OF   130.  YEARS. 
MERCH .   CU.   FT.   -  TREES  6.0   INCHES   D.B.H.    AND  LARGER   TO  4-INCH  TOP. 
BD.   FT.   -   TREES   6.5    INCHES   D.B.H.    AND  LARGER   TO  6-INCH  TOP. 


SAWTIMBER 
VOLUME 
BD.FT. 


-12- 


STANO 
AGE 
(YEARS) 

30 
10 


TREES 
NO. 


1C00 
505 


YIELDS   PER   ACRE   OF   EVEN-AGEO   STANDS  OF   LOOGEPOLE  PINE 

SITE    INDEX  70 
THINNING    I NTFNS I TY-    INITIAL-      120.      SUBSEQUENT-  100. 


ENTIRE   STAND   BEFORE   AND   AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

110 

72 


AVERAGE      AVERAGE  TOTAL 
D.B.H.        HEIGHT  VOLUME 
IN.  FT.  CU.FT. 


4.5 
5.  1 


28 
2<) 


1470 
1010 


"ERCHANT- 
AHLF  VOLUME 
CU.FT . 

210 
210 


SAWT IHBER 
VOLUME 
BD.FT. 


TREES 
NO. 


PERIODIC   INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL  MERCHANT- 
VOLUME     ABLE  VOLUME 


SAWTIMBER 
VOLUME 
BD.FT. 


50 
50 


500 
282 


130 

87 


6.9 
7.5 


41 

42 


2660 
1320 


2000 
1500 


70 
70 


80 
90 


110 
110 


120 
130 


281 
183 


183 
117 


117 

80 


80 
80 


135 
100 


139 
100 

116 

132 
99 

115 

129 


9.4 
10.0 


11.8 
12.5 


14.4 
15.  1 


16.2 
17.2 


55 
56 

61 

65 
66 

69 

72 
73 

75 

73 


3760 
2810 


4470 
3200 


4  580 
3450 


4030 
4720 


3480 
2630 


4220 
3040 


4370 
3300 


3920 
4540 


14500 
11100 


18700 
13700 


20600 
15800 


19200 
22600 


OWARF   MISTLETOE    INFECTION   DID   NOT   OCCUR   DURING   THE   ROTATION  OF   130.  YEARS. 
MERCH.   CU.   FT.   -   TREES   6.0    INCHES   D.B.H.    AND  LARGER   TO  4-INCH  TOP. 
BO.   FT.   -   TREES  6.5    INCHES   C.B.H.    AND  LARGER   TO  6-INCH  TOP. 


STANO 
AGE 
(YEARS) 

30 

30 


50 
50 


70 
70 


TREES 
NO. 


ICOO 
454 


443 
261 


258 
175 


YIELDS   PER   ACRE   OF   EVEN-AGEO  STANDS  OF   LODGEPOLE  PINE 

SITE    INDEX  70 
THINNING    INTENSITY-    INITIAL-      120.      SUBSEQUENT-  80. 


ENTIRE   STAND   BEFORE   AND  AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

110 

36 


62 
84 


107 
77 


AVERAGE 
D.B.H. 
IN. 

4.5 
3.8 


6.3 
6.6 


AVERAGE 
HE IGHT 
FT. 

26 
24 


37 
38 


8.7 
9.0 


51 
51 


TOTAL 
VOLUME 
CU.FT. 

1380 
430 


1760 
1170 


2720 
2000 


MERCHANT- 
ABLE VOLUME 
CU.FT. 

0 
0 


1040 
800 


SAWTIMBER 
VOLUME 
BO. FT. 


2460 
1830 


10100 
7500 


TREES 
NO. 


PERIODIC   INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


90 
90 


110 
110 


120 
130 


172 
113 


111 
77 


76 
73 


116 

80 


110 

80 


11.1 
11.4 


13.5 
13.8 


61 

64 


67 
68 


3510 
2440 


3600 
2620 


3310 
2310 


3430 
2500 


14300 
10100 


15800 
11600 


34 


30 


980 


92  14.9  71  3100  2960  14100 

102  16.0  73  3510  3370  16400 

OWARF   MISTLETOE    INFECTION   STARTED  AT   AGE      10.   AND  RATING  WAS        3.5  AT  AGE  130. 

MERCH.   CU.   FT.   -   TREES   6.0    INCHES  D.B.H.    AND  LARGER   TO  4-INCH  TOP. 

BO.   FT.   -   TREES   6.5    INCHES   D.B.H.  ANO  LARGER   TO  6- INCH  TOP. 


-13- 


STAND 
AGE 
(YEARS) 

30 
30 


TREES 
NO. 


1C00 
454 


YIELDS   PFR   ACRE   OF   EVEN-AGED  STANDS  OF   LODGEPOLE  PINE 

SITE    INDEX  70 
THINNING   INTENSITY-    INITIAL-      120.      SUBSEQUENT-  100. 


ENTIRE   STAND   BEFCRE   AND  AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

1  10 

36 


AVERAGE 
D.B.H. 
IN. 

4.5 
3.8 


AVERAGE 
HE IGHT 
FT. 


TOTAL 
VOLUME 
CU.FT. 

1380 
430 


MERCHANT- 
ABLE VOLUME 
CU.FT. 

0 
0 


SAWTIMBER 
VOLUME 
BD.FT. 

0 

0 


TREES 
NO. 


PERIODIC   INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU.FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWTIMBER 
VOLUME 
BD.FT. 


50 
50 


443 

332 


6.3 
6.5 


1760 
1420 


1040 
930 


70 
70 


329 
235 


127 

95 


8.4 
8.6 


50 
51 


3200 
2420 


2850 
2  180 


11700 

8900 


110 
UO 


120 
130 


227 
227 


198 
198 


177 
153 


134 
134 


156 
156 


10.4 
10.4 


12.0 
12.0 


60 

63 

66 
66 


4020 
4020 


5020 
5020 


3780 
3780 


4750 
4750 


16000 
16C00 


21  100 
21100 


156  12.7  67  5120  4860  22000 

148  13.3  69  4920  4750  21500 

DWARF   MISTLETOE    INFECTION   STARTFO   AT   AGE      10.    AND  RATING  WAS 
NOTE   THAT   NOT   ALL   SCHEDULED   THINMINGS   WERE  POSSIBLE. 
MERCH.   CU.   FT.   -   TREES   6.0    INCHES   D.B.H.    AND  LARGER   TO  4-INCH  TOP 
BD.   FT.   -   TREES   6.5    INCHES   C.B.H.    AND  LARGER   TO  6-INCH  TOP. 


6.0  AT  AGE  130. 


YIELDS   PER   ACRE   OF   EVEN-AGEO   STANDS   OF   LODGEPOLE  PINE 

SITE    INDEX  70 
THINNING    INTENSITY-   INITIAL-      120.      SUBSEQUENT-  80. 

ENTIRE   STAND  BEFORE   AND  AFTER   THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT 

TOTAL 
VOLUME 
.  CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

50 
50 

2500 
2500 

123 
123 

3.0 
3.0 

33 
33 

1830 
1830 

0 
0 

0 
0 

0 

0 

0 

0 

0 

60 

2007 

150 

3.7 

39 

2780 

0 

0 

70 
70 

1565 
1565 

158 
158 

4.3 
4.3 

44 
44 

3360 
3360 

0 
0 

0 
0 

0 

0 

0 

0 

0 

80 

1185 

149 

4.8 

47 

3460 

0 

0 

90 
90 

870 
870 

128 
128 

5.2 
5.2 

50 
50 

3150 
3150 

750 
750 

0 
0 

0 

0 

0 

0 

0 

100 

729 

125 

5.6 

51 

3170 

1  160 

0 

110 
1 10 

611 
611 

120 
120 

6.0 
6.0 

52 
52 

3130 
3130 

1550 
1550 

0 
0 

0 

0 

0 

0 

0 

120 

512 

118 

6.5 

53 

3150 

2060 

0 

130 

429 

115 

7.0 

54 

3120 

2390 

0 

DWARF 

MISTLETOE 

INFECTION 

STARTED 

AT   AGE  10. 

AND  RATING 

WAS 

6.0  AT 

AGE  130. 

NOTE   THAT  NOT  ALL   SCHEDULED   THINNINGS   WERE  POSSIBLE. 

MERCH.  CU.  FT.  -  TREES  6.0  INCHES  D.B.H.  AND  LARGER  TO  4-INCH  TOP. 
BD.   FT.  -  TREES   6.5   INCHES  D.B.H.   AND  LARGER   TO  6-INCH  TOP. 
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STAND 
AGE 
(YEARS) 

50 
50 


TREES 
NO. 


2500 
2500 


YIELDS  PER  ACRE  OF   EVEN-AGEO   STANDS   OF   LOOGEPOLE  PINE 

SITE   INDEX  70 
THINNING   INTENSITY-    INITIAL-      120.      SUBSEQUENT-  100. 


ENTIRE   STANO  BEFORE   AND  AFTER  THINNING 


BASAL 
AREA 
SQ.FT. 

123 
123 


AVERAGE  AVERAGE  TOTAL  MERCHANT- 
D.B.H.        HEIGHT     VOLUME      ABLE  VOLUME 


3.0 
3.0 


33 
33 


1830 
1830 


SAWT I MBER 
VOLUME 
BO. FT. 


TREES 
NO. 


PERIODIC    INTERMEDIATE  CUTS 


BASAL 
AREA 
SQ.FT. 


TOTAL 
VOLUME 
CU. FT. 


MERCHANT- 
ABLE VOLUME 
CU.FT. 


SAWT I MBER 
VOLUME 
BD.FT. 


70 
70 


90 
90 


110 
110 


120 
130 


1565 
1565 


870 
870 


611 
611 


512 
429 


158 
158 


128 
128 


120 
120 


4.3 
4.3 


5.2 
5.2 


6.0 
6.0 


44 
47 


50 
50 


52 
52 


3360 
3360 


3150 
3150 


3130 
3130 


750 
750 


1550 
1550 


118  6.5  53  3150  2060  0 

115  7.0  54  3120  2390  0 

DWARF  MISTLETOE    INFECTION   STARTED  AT  AGE      10.   AND  RATING  WAS 
NOTE  THAT  NOT  ALL   SCHEDULED   THINNINGS   WERE  POSSIBLE. 
MERCH.   CU.   FT.  -   TREES  6.0   INCHES  D.B.H.    AND   LARGER   TO  4-INCH  TOP 
BD.   FT.   -   TREES  6.5   INCHES  D.B.H.   AND  LARGER   TO  6-INCH  TOP. 


6.0   AT   AGE  130. 


Agriculture  —  CSU,  Ft.  Collins 
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